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We study the transmission of botii classical or quantum information through all the phases of a finite XXZ spin 
chain. This characterizes the merit of the different phases in terms of their ability to act as a quantum wire. As 
far as quantum information is concerned, we need only consider the transmission of entanglement as the direct 
transmission of a quantum state is equivalent. The isotropic anti-ferromagnetic spin chain is found to be the 
optimal point of the phase diagram for the transmission of quantum entanglement when one considers both the 
amount of transmitted entanglement, as well as the velocity with which it is transmitted. But this optimal point 
in the phase diagram moves to the Neel phase when decoherence or thermal fluctuations are taken to account. 
This chain may also be able to transfer classical information even when, due to a large magnitude of the noise, 
quantum information is not transmitted at all. For a certain range of anisotropies of the model, a curious feature 
is found in the flow of quantum information inside the chain, namely, a hopping mode of entanglement transfer 
which skips the odd numbered sites. Our predictions will potentially be testable in several physical systems. 



I. INTRODUCTION 

Recently, condensed matter many-body systems have been 
viewed in the light of quantum information. For exarnple, the 
entanglement inherent in them has been investigated ^ . One 
can, however, ask a different question: how does information 
placed on one part of a many -body system pass through such a 
system? Aside its fundamental interest, this question may lead 
to mechanisms for moving information over small distances. 
The idea is to use a finite many-body system such as a spin 
chain (a chain of perpetually interacting stationary spins - a 
one dimensional magnet) as a data-bus |2j]. Many -body dy- 
namics transports information placed on a spin at one end of 
the chain to the spin at its other end with a certain efficiency. 
This is an "all solid-state" bus whose spins and interactions, 
except for those at its very ends, are never controlled. Ap- 
plications could be in moving information between quantum 
registers or for moving classical bits in nano-scale spintron- 
ics. This area, reviewed in jsll, has mainly focussed on per- 
fecting the information transfer (information transmission) by 
clever means - special couplings {4^, encoding |5ij, pulsing 
10], memory 101 etc. 

An interesting question from a condensed matter angle is 
how the above process of information transmission varies with 
the "phase" of the spin chain. By "phase" we mean both the 
form of the spin-spin interactions and the relevant ground state 
resulting from that interaction. In this context, only one study 
has been performed, which involves spin-1 chains fl]. Addi- 
tionally, gapless phases have been shown to be generically bad 
for a "slow" information transmission process that can take 
placej^etween two spins coupled weakly to a many-body sys- 
tem 190. The same slow information transmission process be- 
tween spins coupled weakly to an anti-ferromagnetic (AFM) 
chain has also been studied ifioll . However, there is no inves- 
tigation yet of information transmission as a function of the 
phases of the simplest, namely the spin- 1/2 chain, when all 
spins are coupled equally strongly so that information trans- 
mission is fast. Instead, a majority of the work has simply 
assumed a fully polarized (symmetry broken) ferromagnetic 
(FM) initial state of the spin chain Here we study the pro- 



cess of information transfer through all phases of a S* = 1/2 
XXZ Heisenberg-Ising chain which models a range of realis- 
tic materials and, according to Ref. jU], is the most impor- 
tant paradigm in low-dimensional quantum magnetism. Us- 
ing finite chains (the case relevant for information transmis- 
sion) and exact diagonalization, we identify the point in the 
phase diagram which provides the optimal data-bus in absence 
of any encoding, engineering, control etc. Interestingly, this 
turns out to be the "isotropic" AFM phase which is the most 
interesting phase ifTlll of the XXZ model. Here the ground 
state has complete SU(2) symmetry and contains significant 
"quantum" correlations or entanglement. This phase is, per- 
haps, also the most common, as it appears in the ubiquitous 
Hubbard model at strong repulsion and half filling. Addition- 
ally most solid-state spin chains such as the famous KCuFa 
ifllll . engineered atomic-scale spin chains jl2j] and doped ful- 
lerine Sc@C82 chains iH] are naturally AFM. 

This study is an example of non-equilibrium dynamics in 
many -body systems, currently a topic of intense activity ifTill . 
Our dynamics is induced by suddenly coupling a single spin 
(the one bearing the information) with one end of a finite spin 
chain. For a range of phases, certain spin correlation func- 
tions behave curiously during this dynamics so that the initial 
state of the added spin hops through the chain skipping al- 
ternate sites. Additionally, information transmission exhibits 
contrasting behavior in the FM and AFM parts of the so-called 
XY phase and has a sharp jump at the boundary of the XY and 
FM phases. 

The structure of this paper is as follows: in section 
we show that entanglement distribution through an arbitrary 
channel is equivalent to the process of transferring a quan- 
tum state through the channel. In section |III| we introduce our 
model, i.e., XXZ Hamiltonian, and in sectionHVlwe consider 
the entanglement distribution via whole phase diagram of a 
XXZ chain. This is followed by an explanation in sectionlV] 
In section[Vl]we characterize the effect of the channel. In sec- 
tion lVIII and section lVlIIl the thermal fluctuations and interac- 
tions with a bath are investigated respectively. Classical com- 
munication through this system is the subject of section HXl 
which is followed by considering the information flow "inside 
the chain" in section|X] In section IXH we give some potential 
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physical realizations which might test our results, while we 
summarize our results in section lXlIl 



n. EQUIVALENCE OF STATE TRANSFERRING AND 
TELEPORTATION MODELS OF INFORMATION 
TRANSMISSION 

In order to transfer information from one place to another 
we have to transfer a state (say the state of a spin) which en- 
codes some information. In particular, when we are thinking 
about quantum information transmission, to quantify the qual- 
ity of transmission, we compute the fidelity between the sent 
and the received state. Since this fidelity is dependent on the 
initial state it is preferable to take the average value of the 
fidelity over all possible equi-probable initial states. This av- 
erage fidelity makes it possible to compare transmission qual- 
ity of different channels and different schemes of information 
transmission. To send quantum information from sender to 
receiver one can think about two different strategies. In the 
first strategy, which is called "quantum state transferring", the 
quantum state is sent through the channel directly. Because 
of the interaction between the channel and quantum state they 
become entangled and state transferring is imperfect in the 
sense that the fidelity between the received state and the ini- 
tial state is less than one. On the other hand, instead of using 
state transferring one can use teleportation for sending quan- 
tum information. Teleportation is based on a shared entan- 
gled pair between sender and receiver which plays the role of 
the resource ifisll . In this second strategy of sending quan- 
tum information, the sender generates a maximally entangled 
pair, keeps one part, and send the other one to the receiver 
through the channel. This shares an entangled pair between 
both sides of the channel and teleportation between sender and 
receiver can be used for information transmission. However, 
this fact that the entanglement of the shared pair is not maxi- 
mal makes the teleportation imperfect. The importance of the 
second strategy is that we just send one part of the singlet state 
through the channel and it is not necessary to study the effect 
of the channel on an arbitrary state. What we show in this sec- 
tion is that the average fidelity in both strategies are the same. 
This was already shown in llql using a different technique for 
arbitrary dimensions of the Hilbert spaces and here we prove 
it again, just for qubits, using a much simpler language. 

Let's start with the state transferring. In this case quantum 
state goes through the channel. An arbitrary quantum channel 
^ is completely determined by a set of Kraus operators {K,n} 
such that the output of the channel is 

m m 

where ps is the input state of the channel, p,. is the output 
state received by the receiver and ST stands for "State Trans- 
ferring". Here we start from the most general form of a qubit 
state lips) = cos6'/2|0) + e*"^ siii6l/2|l) as the input. Af- 
ter interacting the pure input state ps = \tps){i^s\ with the 
channel the output state (given by Eq. ([T])) is generally a 
mixed state. Fidelity between the received and the sent state 



is easily computed as F^'^{0,(j)) = (-^slpf^lV'.?) which is 
dependent on input parameters and 0. To get an input in- 
dependent quantity we average the fidelity over all possible 
input states, i.e. the surface of the Bloch sphere, with uniform 
weight. With a straight forward computation we end up with 

= ^ / F'^{0,4>)fi\nededct> = i + i^ |Tr(if„)p, 

m 

(2) 

where Tr{.) = Trace{.). 

Now, we try to use the teleportation strategy for sending 
quantum information. To achieve this strategy we prepare a 
pair of singlet state 

I*-) ^ M> . (3) 

Then we keep one part of the pair in the sender and send the 
other part through the channel ^. Since the first part in the 
sender does not interact with the channel the whole effect of 
the channel is explained by 

Pout ^i® m~){i^~ I) = E ^ ® ^"'1'/'") 1^ ® 

m 

(4) 

Generally, the output state pout is not a maximally entangled 
state so when it is used as the resource of the standard tele- 
portation scheme ifTsIl it gives an imperfect teleportation in 
the sense that the final achievable fidelity is less than one. In 
IitIi it has been shown that teleporting p^ with using noisy re- 
source Pout generates the following state as the output of the 
teleportation. 

3 

P^^ = ^ Tr{poutEm)(^raPs(^m, (5) 
m=0 

where TP stands for "Teleportation", Em = crm\'^^){'4^ — 
\am and (7^ are Pauli matrices (ctq = I,<^i.2,3 = <^x,y.z)- 
Similar to the first strategy, fidelity of the received and the sent 
state is defined as F'^^{9,4i) = {ips\p'^^H's) and average 
fidelity for input states is easily computed over the surface of 
the Bloch sphere. The average fidelity of teleportation scheme 
is 

Faf = ^/ F^^(0,0)sin(0)d0d^ 
1 ^ 

= Tr{Eopout) + 3 X! FriEmPout) 

ni—1 

1 + 2Tr{EoPout) 
= ^ . (6) 

The parameter Tr(i?oPoiit) = {4'~\Pout\'>P~) is called singlet 
fraction and as it is clear from Eq. (|6]l that it completely cap- 
tures the quality of the transmission. It is also clear from Eq. 
^ that to have an average fidelity above 2/3, which is acces- 
sible to the classical teleportation, singlet fraction should ex- 
ceed 1/2. Using the form of pout in Eq. ^ and expanding the 
singlet as Eq. @ one gets Tr{EoPout) = j Em \Tr{Km)\^. 
Substituting this value in Eq. Q shows that = F^J . 
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Getting identical average fidelity in both strategies is a very 
important result in quantum communication which shows the 
average effect of a channel can be captured just by transferring 
one part of the singlet state through the channel and comput- 
ing the singlet fraction. However, sharing an entangled pair 
between sender and receiver has an advantage, namely that 
after a few transmissions the total (generally noisy) entangle- 
ment can be converted by local actions [18] to nearly a pure 
singlet. This can be used to transmit any state near perfectly 
using quantum teleportation. So, because of the importance 
of the amount of entanglement shared between the the sender 
and receiver, and its above proven equivalence to the more 
straightforward transmission of quantum states, we mainly fo- 
cus on the entanglement distribution through the phase dia- 
gram of the XX Z Hamiltonian. 



III. INTRODUCING THE MODEL 

We consider a spin chain as a channel for information trans- 
ferring and we study the property of each phase of the chain 
on the quality of information transmission. We take one of the 
most well known models in condensed matter physics, namely 
XX Z spin chain. The Hamiltonian of the open XXZ chain of 
length Nch is 

H,,^J J2 {«+i + + Aafaf+ J, (7) 

with J being a coupling constant, A being the anisotropy and 
fj^-y^^ being Pauli matrices for site k. This Hamiltonian has a 
rich phase diagram. For A = 1 and J < this interaction is 
the FM Heisenberg chain widely discussed in the context of 
quantum communication 10, S Hi ■ More interesting regimes 
exist for J > and different values of A UH]. A < -1 
is the FM phase with a simple separable biased ground state 
with all spins aligned to the same direction. — 1 < A < 1 
is called XY phase, which is a gapless phase and consist of 
two different legs, ferromagnetic half (— 1 < A < 0) and anti 
ferromagnetic part (0 < A < 1). 1 < A is called Neel phase, 
where the spectrum is gapped and we get nonzero staggered 
magnetization. In the limit A ^ 1 it takes the form of Neel 
states (1010101. ..01)). 



IV. INFORMATION TRANSMISSION THROUGH WHOLE 
PHASE DIAGRAM OF XXZ HAMILTONIAN 

Though information transmission can be investigated either 
classically or quantum mechanically we will primarily exam- 
ine quantum information transmission and devote one section 
later to classical information transmission in the same sys- 
tems. Of course, the most natural setting would be sending 
the state of a single spin through the chain. However, be- 
cause of previous "equivalency" discussions in section (Hil l, 
we will examine the transmission of one part of a two spin 
maximally entangled state of the form (|3]l while a spin chain 



channel (spins 1 to Nch) is in its ground state of some Hamil- 
tonian Hch- At time t = 0, the interaction of the 0th and the 
1st spin is suddenly switched on while 0' is kept isolated from 
the rest. The ensuing dynamics transports the initial state of 
the 0th spin through the chain to the -/Vc^th spin with some 
efficiency, so that after a while 0' will be entangled with Nch- 
As the singlet has the same representation in any bases, the 
above entanglement transfer already subsumes within it "state 
transfer in arbitrary basis" and is thus very general. 

The reader may naturally question how general the above 
physical setting (couplings etc.) of transferring entanglement 
through a spin chain channel is. Indeed one could have taken 
weaker or stronger or different couplings at the sending and 
receiving ends. However, weaker couplings generally lead to 
"slow" transfer schemes which will be susceptible to deco- 
herence. On the other hand, if we really do have stronger 
couplings or different couplings available at our disposal, we 
could just use them for the whole chain for faster and poten- 
tially better transfer, rather than using those special couplings 
only at the ends. So we think that the most natural question 
to investigate is to simply place a spin encoding the unknown 
state to be transmitted at one end of the chain, and couple 
it with the same coupling as present in the rest of the chain 
(which, as we know from the previous section, is equivalent 
to the type of entanglement transmission considered by us). 
In any case, without putting some restrictions on the coupling 
model at the ends, there is too much freedom in the prob- 
lem, and it may not be possible to give a precise answer to 
the effectiveness of a phase to transfer quantum information. 
Moreover, also note that we are not considering the genera- 
tion of entanglement from inside the spin chain, which is an 
altogether different problem lfl9ll . but merely the transmission 
of entanglement through the chain. 

Note that for our scheme we require the chain initially in a 
unique ground state \ipg)ch and this may have to be selected 
out by applying an arbitrarily small magnetic field (for odd 
Nch AFM chain and the FM chain). The interaction between 
the 0th and the 1st spins (the interaction turned on at t = 0) of 
the channel is assumed to be of the same form and strength as 
the rest of the interactions, namely 

Hi = J{a^af + a^af + Aa'^'^!), (8) 

With the O'O singlet, the total length of the system consid- 
ered is thus N = Nch + 2 with the initial state being 

|i^(0)) = |V")o'o® |^s)ch, (9) 
and the total Hamiltonian being 

H = h,® {Hch + Hi). (10) 

So that 0' never interacts with the rest. Also note that H is 
simply a Hamiltonian of a single spin chain O...Nch of length 
+ 1. As the aim is entanglement distribution, we are in- 
terested at the times that the entanglement between spins 0' 
and Nch peaks. By turning on the interaction between spin 
and spin 1 of the channel the initial state evolves to the 
state iV'(i)) = e^'^*|?/'(0)) and one can compute the den- 
sity matrix pij ~ t'>'ij{\ip(t)){'4-'{'t)\} where the meaning of 
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trfj is the trace over whole of the system except sites i and 
j (We fix i = 0' in this paper). The general form of a two 
spin density matrix pij in X X Z systems in the computational 

(joO), |01), 1 10), 111)) basis are 
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where all the elements of the matrix are real and they can be 
written in terms of one and two point correlations, 
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(af(t)a|(t))T(af(t))±K(i))} 



(12) 



where cr"(i) = e'^V^e^'^* (for all a = x,y,z) is the 
Heisenberg picture of cr" and <> means expectation value 
according to the initial state (|9]l. The concurrence as a mea- 
sure of entanglement ll2lll for this general density matrix (fTTT i 
is £^ = 2max(0, \z\ — yu+ir^ ) which is a function of time 
dependent correlators and expectation values. Noteworthy, is 
that for A > — 1, when the initial state of the channel is not 
symmetry broken, then symmetry considerations and the fact 
that and 0' are initially anti-correlated in a singlet, imply that 
the entanglement between 0' and j can be written as 



Eo>, = max(0, |(ai?(0)aj(t))|--(a5(0)a|(i)) 



(13) 



which is solely written in terms of the two-time correlation 
functions of the spin chain O...Nch- It should be noticed that 
though two point correlations of the XX Z Hamiltonian have 
been studied intensively in the literature and their asymptotic 
behavior is known, the correlations here are different since 
they are computed in terms of the initial state (|9]l which is 
not the ground state of H. In addition, if one ignores (traces 
out) spin 0', our study can be regarded as an analysis of two 
time correlation functions during the non-equilibrium dynam- 
ics that ensues when the interaction of a spin in a random state 
with one end of a spin chain is switched on. Singlet fraction 
of the state po' n, which was shown that is directly related to 
the average fidelity of state transferring, can be computed also 
from pqi TV easily 
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In Fig. [T]we have plotted both entanglement and singlet 
fraction of po'iv in terms of time for a particular point in the 
phase diagram, namely A = 1. As it can be seen from fig- 
ure, singlet fraction always oscillates while the entanglement 
just peaks at certain times which we call optimal time topt- 
When entanglement peaks, singlet fraction also has a peak 
which shows that final state is more similar to singlet than 
other Bell states. The time that one can afford to wait for the 
entanglement between 0' and Nch to attain a peak is restricted 
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FIG. 1 : (Color online) Entanglement and singlet fraction in terms of 
time for a chain of length A'' = 20 and A = 1. 
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FIG. 2: ((Color online) a) Attainable entanglement in the first peak in 
terms of A for different lengths (J=l). Inset shows the optimal time 
topt that the peak happens during the evolution, b) Singlet fraction 
F at topt in whole phase diagram. 



by practical considerations such as the decoherence time, re- 
quired speed of connections in a quantum network etc. So we 
restrict ourselves to the first peak of the entanglement in time. 
To compare the performance of different phases of the Hamil- 
tonian (|7]i in transferring the entanglement, we have plotted 
the amount of entanglement in its first peak {t = topt) in terms 
of anisotropy A for the chains with different lengths in Fig. 
and the associated singlet fraction at the same time in Fig. ^p. 
One interesting feature is that the entanglement transmitted 
dips on the XY side of A = — 1 and sharply rises on its FM 
side which captures the first order phase transition at this point 
(that such a change in behavior is seen despite the finite size 
is interesting). In addition, this transition is also marked by 
a steep rise in the time required to reach the first peak in en- 
tanglement. This is shown in the inset to Fig. |2^, which also 
shows that the speed at which entanglement is propagated in- 
creases monotonically in the XY regime as one goes from 
A = —1 to 1. This fact is commensurate with the spin wave 
velocity increasing with A as sin(cos~^ A)/ cos~^ A in this 
regime jH] which we will discuss it in more detail later. In the 
XY phase we can recognize two distinct regimes. In its FM 
sector (— 1 < A < 0) entanglement falls rapidly by decreas- 
ing A, while in its AFM sector (0 < A < 1) the entanglement 
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FIG. 3: (Color online) a) Optimal time topt for both FM (J — — 1) 
and AFM (J = +1) in terms of length A^. b) Entanglement at the 
first peak versus the length A'' for both FM (J = —1) and AFM 
(J = +1) chains when A = 1. 



is always good and increases by increasing A. After A = 1, 
when the transition from the XY to the Neel phase happens, 
the entanglement starts falling with increasing A, as the Ising 
term (t|(t|^j^ dominates which, by itself, does not transfer en- 
tanglement. Note also a subtle feature that even outside the 
XY regime, for A > 1, entanglement falls much slower with 
|A| than for A < —1. In general, AFMs are thus better, even 
with similar degrees of anisotropy. Fig. [2h also shows that the 
isotropic AFM Heisenberg interaction (A = 1) not only is the 
best for transferring the highest amount of entanglement in the 
entire phase diagram, but also it has the highest speed in the 
XY phase. In Fig. I2J), where singlet fraction F has been plot- 
ted in the entire of phase diagram, in the FM phase (A < — 1) 
singlet fraction is always less than 1/2 which shows despite 
the fact that entanglement is non-zero, quantum communica- 
tion has no benefit over classical communication. Same thing 
will happen in the Neel phase when F becomes less than 1/2 
for quite large A's. 

The effect of the length of the chain on the quality of trans- 
mission has been shown in Fig. |3] We only concentrate on the 
best A = 1 (isotropic AFM) point, as it is the best point in the 
phase diagram, and we compare the results with FM chains 
which have been predominantly studied so far In Fig. [3^ we 
have plotted the time at which the first peak in entanglement 
for different lengths. It is clear that the speed of entangle- 
ment transmission through the AFM (J > 0) chain is higher 
than FM (J < 0) chain independent of the length. In Fig. [Sji 
the amount of entanglement in the first peak has been com- 
pared for both the AFM and FM cases, from which it is clear 
the the entanglement transmitted in the case of AFM chain 
has a distinctively higher value irrespective of length. Note 
also a visible even-odd effect on the amount of entanglement 
transmitted (and hence on two-time correlations) which will 
be interesting to observe in finite chains. 

We have shown that in absence of any of the sophisticated 
techniques for perfecting spin chain communications, which 
come at a price, and may be hard to implement and also if 
one wanted to transfer information fast (i.e, refrain from very 
weak couplings), the isotropic AFM is the best channel in 



the entire phase diagram of the XXZ chain. We now esti- 
mate the efficiency with which local processing at the oppo- 
site ends of the spin chain and classical communication be- 
tween them (a process called entanglement distillation fWi) 
can establish a nearly perfect singlet for an isotropic AFM 
channel. For instance by using the recurrence algorithm for 
distillation ifisll in a chain of length 10, for which entangle- 
ment E = 0.8638, starting from 9 impure pairs on average 
leads to a nearly singlet state with entanglement E ~ 0.9920 
after 7 iterations, and for a chain of length 20, for which entan- 
glement E = 0.7162, we need to start with 17 impure pairs 
to get a singlet state with entanglement E = 0.9926 after 
9 iterations. This perfect singlet can then be used for send- 
ing quantum states perfectly through teleportation. It is worth 
pointing out that in different phases the spin chains represent 
different types of quantum channels. While in the FM phase, 
it is known to be an amplitude damping channel (transmits 
|0) and |1) asymmetrically 1I2II), the A = 1 point affects a so 
called depolarizing channel (also noted in lUOll ') where po' n^h 
is the mixture of the singlet state and the Identity. Such po' n^^. 
is particularly suited for distillation protocols llSll . 



V. EXPLANATION 

When the phase of the system changes, not only the Hamil- 
tonian causing the time evolution varies, but also the ground 
state, and consequently the initial state (|9]l varies, and we have 
a different behavior for information transmission through the 
chain. Results of the Fig. [2^ shows a dramatic and discontin- 
uous change of entanglement at A = —1 which is related to 
a first order phase transition at this point and two completely 
different class of ground states of the XY and the ferromag- 
netic phase. At point A = +1 entanglement falls continu- 
ously when we go from XY phase to the Neel phase. This 
continuous change represents a second order phase transition 
at this point. Beside these two phase transition points there is 
a sharp drop of entanglement around A = —0.5 which is very 
peculiar since there is no phase transition at this point. Also 
from the inset of Fig. |2^ it is clear that the optimal time which 
one has to wait to get a peak goes up drastically for A < —0.5. 
This strange property inside the XY phase is certainly not be- 
cause of a phase transition. The reason of this slow dynamics 
and bad transmission is hidden behind an intrinsic property of 
the spin chain namely, the spin wave velocity. 

Field theoretic techniques have been used to capture the 
asymptotic behavior of correlation functions in spin chains 
12311 . For the general XXZ Hamiltonian, it fails to get all 
prefactors and exact solutions, but it is able to get the qual- 
itative behavior of correlations successfully in the thermody- 
namic limit. Correlation functions in our problem are different 
from those obtained by field theory in at least two ways. First 
of all, we consider very finite chains, since the idea of using 
spin chains as quantum channels is valuable only for a finite 
distance. Secondly, all the dynamical correlations which are 
computed asymptotically are associated with the ground state 
of the system while in our problem correlations are computed 
for the initial state (|9]l which is not the ground state. Despite 
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FIG. 4: (Color online) 1/vf (for infinite chain) and topt (for a chain 
of length iV = 20) versus A. 

these differences, we still can use some well known results 
of the field theoretic techniques. For example, the dynami- 
cal correlation functions in the XY phase (—1 < A < 1) 
in the asymptotic thermodynamical limit have the following 
form [21; 

K(0K(t))^ (-1)1^-^1 (15) 

where a = x,y,z and 

= Vy = ^iVz = I ■ (16) 

TT 

Moreover, in Eq. (fTST i is the spin-wave velocity (this 
quantifies the propagation velocity of excitations in the chain) 
which has the following form 

sin(cos^^A) 

VF(X — ^ (17) 

cos~^ A 

Unfortunately, the above asymptotic forms of ((T"(0)(T^(t)), 
validfor |j — fc| >> are i/ngM/ar specifically at |j — /ej ^ 
vpt, which is the regime relevant to optimal quantum commu- 
nication from jth to the fcth site (i.e., when the information, 
possibly traveling at a velocity vp reaches its destination). So 
one can only use some aspects reliably from the above formu- 
lae. One of this is the velocity vp of propagation of the cor- 
relations (and hence information). In Fig. |4]we have plotted 
\/vf (which is for an infinite chain) and topt (for a chain of 
length N = 20) in terms of A. As Fig. |4]clearly shows, both 
of these quantities behave in a strikingly similar way (the gap 
between two curves is not important since one can multiply 
them by some constants). As a consequence, for A < —0.5 
the propagation velocity is very slow, and one has to wait a 
long time to receive some information at the other side of the 
chain. This very slow dynamics also means that we will get 
a sharp fall in the entanglement as well as all other quantities 
which propagate through the chain unless we are willing to 
wait for very very long times. 

As far as the question of why the isotropic (A = 1 point) 
is the best point in the phase diagram in terms of a maximum 
of entanglement, recall that entanglement is given in terms of 



some dynamical correlation functions as in Eq. (fTsT l. Remem- 
ber though, that these correlations are not evaluated for the 
ground state, so cannot be strictly substituted by the known 
dynamical correlation functions to get any quantitative infor- 
mation, but perhaps only a qualitative picture, as we discuss. 
One can see from Eq.lfTSt that the A = 1 point is the best for 
the propagation of correlations along the z direction. At points 
with A < 1 the z component of correlations does not propa- 
gate as well as the x component, and, in fact, near to A = — 1 
it is expected not to propagate at all (t]2 « oo). Thus the term 
— ^{(jQ{0)aj{t)) in Eq. ( fT3] l for entanglement, which is pos- 
itive, contributes more and more as we approach A = 1 and 
gives a higher entanglement. It is true that as we approach the 
isotropic point from A < 1 side, the r]x rises (i.e., propaga- 
tion of correlations in the x direction deteriorates somewhat). 
However, it must be that the gain from the better propagation 
of correlations in the z direction more than compensates for 
the deterioration of the propagation of correlations along the 
X direction. The reason is that rj^ changes from oo to 1 (huge 
gain), while rj^ only goes from to 1. 

As far as the intriguing dip after A = —0.5 is concerned, 
we are not yet in a position to explain it. It seems that the 
behavior expected at A —1 where both the velocity of 
correlations and their propagation quality along the z direction 
are worst, starts to happen quite a bit before the actual point. 



VI. CHANNEL CHARACTERIZATION AND EVEN-ODD 
EFFECT 

As it is clear from Fig. in the case of AFM chain, en- 
tanglement has a zig-zag behavior when N varies while it 
behaves uniformly for FM chains. This even-odd effect for 
AFM chains has a fundamental reason. In even chains when 
A > — 1 the total magnetization of the ground state is always 
zero and because of the rotational symmetry in the ground 
state one can exchange all |0)'s and |l)'s while the ground 
state remains unchanged. In other words, in even chains for 
A > — 1 we always have 

Crf''^'^^g)ch = \^g)ch. (18) 

This symmetry, which is absent in FM chains and also in 
each of the doubly degenerate ground states of the odd 
chains, has a profound effect on the transmission char- 
acteristics of the chain. In even chains, the effect of 
the chain is completely recognized by the Kraus operators 
{\/PiI,^/P^-'^x,^fp^'7y,^/pI<7z}- Where, pi^x,y,z are posi- 
tive and their summation is equal to 1 so, one can explain 
the effect of this channel such that it applies one of the Pauli 
operators (including identity) with some probability to the in- 
put state. Obviously, these probabilities are dependent on the 
length A^, time t and anisotropy A. So using the above Kraus 
operators we get the following form for the state po' N^h. 

Po'N^, = Pi{tM-){4r\+p^{t)\q^-){(^-\ 

+ (19) 
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FIG. 5: (Color online) a) Entanglement in terms of temperature in 
a chain of length A'^ = 10 for the isotropic case (A = 1) in both 
FM (J = -1) and AFM (J = +1) phase, b) Entanglement in 
whole phase diagram for different temperatures in the chain of length 
N = 10. 



where, 



1^^ 



±, _ |01)±|10) 



V2 

±, _ |oo)±|ii) 



^/2 



(20) 



are Bell states. Thus it means that po'Nah is diagonalaized 
in the Bell basis. This channel is called Pauli channel in the 
literature. Since in the Hamiltonian ( fTOl i there is no differ- 
ence between x and y directions we always have p^, = py for 
XX Z chain. At the point A = 1 where all directions become 
identical we have Px = Py ^ Pz and channel is the famous 
depolarizing channel. 

For the case of odd N, characterization of the channel is not 
yet known and one can just consider it numerically. Since in 
odd chains the ground state of the system is degenerate, to take 
one of them we apply a small magnetic field in the z direction 
to break the symmetry. In this case the total magnetization is 
±1 (dependent on the direction of the magnetic field: ±z) and 
the symmetry ( fTST i does not hold anymore. 

For A < — 1 since the ground state is ferromagnetic and all 
spins are aligned the type of the channel is amplitude damping 
II2II so then the even-odd effect vanishes and channel behaves 
uniformly for all N. It worths to mention that in entangle- 
ment distillation procedures ifisll Werner states, which are a 
mixture of Bell states, are distilled more easily than the other 
states 1 180 so transferring the singlets through a Pauli channel 
has this advantage that the final state is very close to a Werner 
state (at A = 1 it is exactly a Werner state) and one can dis- 
till them more easily than those which are gained through the 
transmission of other channels such as amplitude damping. 



VII. THERMAL FLUCTUATIONS 

Generally, when system is in non zero temperature, the state 
of the channel before evolution is described by a thermal state 
- — instead of the ground state, where (3 = 1/KbT and 



FIG. 6: (Color online) a) Entanglement in terms 7 in a chain of length 
= 8 for the isotropic case (A = 1) in both FM (J = -1) and 
AFM (J = +1) phase, b) Entanglement in whole phase diagram for 
different noise strength 7 in the chain of length N = 8. 



Z is the partition function. So in this case the initial state of 
the system is 



p(0) = |Vj-)(^-|® 



(21) 



We assume that the thermalization time-scale of the system 
is large so that one can consider the unitary dynamics start- 
ing the initial state ( I2TI ). So, after time t system evolves to 
p{t) = Up{0)U^ and the target state po'Wch (0 '^^n be gained 
as before by tracing out the bulk of the chain po' n^h (^) = 
tr^,^ ^ {/)(<)}. Entanglement of the state po' n^^ (^) opti- 
mal time has been plotted in Fig. in terms of initial temper- 
ature for both FM ( J = -1) and AFM chains ( J = +1). As 
it is clear from the figure, increasing the temperature always 
destroys the entanglement but it has less effect on AFM chain. 

In Fig. the entanglement in whole phase diagram has 
been plotted for different temperatures. When temperature 
rises entanglement survives more for fully symmetric Heisen- 
berg point (A — 1). Specially FM phase is highly sensitive 
to thermal fluctuations, and entanglement is destroyed rapidly 
when temperature rises. Furthermore, we found that, optimal 
time topt which the entanglement peaks is almost independent 
of the temperature and varies very slowly in the entire of phase 
diagram. 



VIII. INTERACTION WITH BATH AND DECOHERENCE 
EFFECT 

In practical situations it is impossible to isolate a quantum 
system from its environment. In the case of Markovian inter- 
action between the system and the environment, a Lindblad 
equation describes the evolution of the system 



(22) 



where £{p) is the Markovian evolution of the state p. Let us as- 
sume an environment which has no preferred direction. Even- 
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tually the interaction should have the following form 

^(p) = -|EE{/'-'^"'°^"}' (23) 

i a. 

where index i takes 0', 0, ...,Nch and a gets x, y, z, and the 
coefficient 7 stands for the rate of decoherence. In Fig. [6^ 
we have plotted the entanglement in terms of noise strength 
7 for both FM (J = -1) and AFM (J = +1) chain. Since 
the figure clearly shows entanglement decays exponentially 
by increasing 7 but like the thermal effect AFM chain is more 
resistive against Markovian noise. In Fig. I6J), we have plotted 
the attained entanglement in when A varies in whole phase 
diagram. As the figure shows the noise effect always kills the 
entanglement and similar to the thermal fluctuations optimal 
time is almost independent of noise parameter 7. Surprisingly 
the best point in the phase diagram is not the case of A = 1 
and it moves down inside the Neel phase. The reason of this 
interesting phenomena comes from the velocity of dynamics. 
Since the evolution in the Neel phase is faster (see the inset of 
Fig. |2^) and entanglement peaks earlier, decoherence has less 
opportunity to interfere and shows its destructive effect. 

IX. CLASSICAL COMMUNICATION 

We now comment on the classical information transmission 
through spin chains, which may be interesting for spintron- 
ics. To quantify the amount of classical information which 
each channel can transmit concept of the classical capacity 
has been introduced. The classical capacity of the channel ^ 
(introduced by the Kraus operators ([T]) in a very general case) 
gives the maximum amount of classical information that can 
be reliably transmitted per channel use. In calculating the 
classical capacity it is necessary to perform a maximization 
over multiple uses of the channel 

C = maXn——, (24) 
n 

where C„ is the classical capacity of the channel ^ which can 
be achieved if the sender is allowed to encode the informa- 
tion on codewords which are entangled only up to n-parallel 
channel uses. The value of C„ is obtained by maximizing the 
Holevo information I24<1 at the output of n parallel channel 
uses, over all possible input ensembles {pi, pi}, i.e. 

Cn = max{p.^p^}Hn{^'^",{pt, Pi}), (25) 

where Hn{£,®", {pi, Pi}) is the Holevo information which is 
defined as 

Hn = {^(^®"(Ep^A.)) -Y.P^e''^P^)}■ (26) 
i i 

Here p,;'s are probabilities, p^'s are n-qubit codewords (ei- 
ther entangled or separable) and S is the von Neumann en- 
tropy. Unfortunately computing the classical capacity is an 
extremely hard task since it needs a very difficult maximiza- 
tion. But recently the classical capacity of the depolarizing 



channel has been computed ||25ll and it was shown that this ca- 
pacity can be achieved by encoding messages as products of 
pure states belonging to an orthogonal basis, and using mea- 
surements which are products of projections onto this same 
orthogonal basis. So due to the fact that entanglement does 
not increase the capacity of the depolarizing channel all maxi- 
mization shrinks to compute the single shot capacity Ci as the 
real capacity of the channel. 

As we discussed in section (lIVI i. for even chains XXZ 
Hamiltonian is a Pauli channel. At isotropic point (A = 1) 
it is a depolarizing channel, which Ci is the real capacity and 
entangled inputs do not increase it. This motivates us to study 
single-shot classical capacity of the XXZ Hamiltonian for 
pure orthogonal input states. However the single-shot capac- 
ity which is computed over pure orthogonal input states is not 
necessarily the real capacity of the channel (except at the point 
A = 1) but at least it gives us a lower bound of the classical 
capacity. To have the form of Pauli channel, we also restrict 
our study just to the even chains. 

We start with the most general form of the orthogonal pure 
qubit states 

IV'i) = cos^|0)+e'^sin^|l) 

IV^2) = sin^|0)-e'^cos^|l), (27) 

where {) < 9 < tt and < (/) < 27r. For input ensemble, we 
associate the probability pi to the input state pi ~ 
and similarly probability p2 to the state p2 — \4'2) {4'2\- When 
each of these states goes through the channel we get 

i{Pi) = Pi Pt + PxO-xPiO-x + PyO-yPiCTy + PzazPidz, I 1,2. 

(28) 

Where in the above equation pi^x,y,z are dependent on time t 
and anisotropy A. It is easy to see that S'(C(pi)) and S{£^{p2)) 
are equal and independent of (j>. Thus the second term in 
the Holevo information ( |26] ) is PiS{£,{pi)) + P2<S'(C(P2)) = 
S{£,{pi)), which is independent of pi and p2 and it is just 
dependent on 9. We can easily maximize the first term in 
the Holevo information ( |26] l for all values of 9 by choosing 
Pi = P2 = 1/2 such that S{^{pipi +P2P2)) = 1 so to maxi- 
mize the Holevo information Hi one should just find 9 = 9 opt 
such that minimize S{£,{pi)). 

Our anaytic computation shows that 

«/ Pz > Px ■ Oopt = or TT. 

if Pz <Px ■ Oopt = 7r/2. 

if Pz =Px ■ Oopt = arbitrary, (29) 

where the situation pz = Px is associated to the depolarizing 
channel (A = 1) which for any value of < 6* < tt the clas- 
sical capacity is achieved. Important point is that the optimal 
input ensemble is independent of phase (p which gives us a lot 
of degrees of freedom for input states. This also was expected 
due the symmetry of the x and y directions in our Hamiltonian 

In Fig. |7^, we have plotted the classical capacity in terms 
of A. This figure clearly shows that Hi is quite flat in the XY 
phase and it suddenly falls for A < —0.5. More interesting 
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FIG. 7: (Color online) a) Holevo information Hi in terms of A for 
different lengths, b) Optimal angel 6opt for input set of states in terms 
of A. 
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FIG. 8: (Color online) a) Entanglement between site 0' and other 
sites in the chain during the evolution in an AFM chain (A = 1) 
of length = 6. b) One configuration of the states which are not 
accessible energetically. 



result has been shown in Fig. [TJ? where optimal 9 has been 
plotted in whole phase diagram. It shows that when we cross 
the point A = 1 from XY phase to the Neel phase suddenly 
the optimal ensemble changes from orthogonal states on the 
equator {6 = 7r/2) of the Bloch sphere to the states on the 
poles (9 = 0). Within the XY phase, around A = -0.35 the 
optimal input ensemble changes such that optimal states for 
-1< A < -0.35 are gained by 61 = and for -0.35 < A < 
1 are obtained by = tt/ 2. In the the FM phase (A < -1) 
the transmission is completely different and it is explained as 
an amplitude damping channel j^. For this channel it was 
shown that Ci is achieved by the inputs given as Eq. ( l27b for 
6* = 7r/2 ||26ll. 

It is interesting to check the classical capacity of the chan- 
nel when it can not transmit quantum information. So for a 
chain of length = 8 with the noise parameter 7 = 0.3, en- 
tanglement can not be transferred because of the large noise 
(quantum information transmission is impossible) but at opti- 
mal times for isotropic case (A = 1) one gains Ci = 0.3931 
for AFM chain (J ^ +1) and Ci 0.1453 for FM chain 
(J = -l). 



X. ENTANGLEMENT PROPAGATION THROUGH THE 
CHAIN 

A curious feature emerges in the propagation of entangle- 
ment through chains with even numbers of spins. For A > 0, 
there is never any entanglement at any time between site 0' 
and odd sites and entanglement seems to hop through the 
chain. If one takes an approach whereby one draws a bond 
for the presence of strong entanglement and dotted for very 
weak entanglement (< 0.1), the open ended ground state will 
be depicted as a dimer (remember it is not an exact dimer) 
iIztIi . Appending a singlet of spins and 0' at one end of the 
chain, makes the total system look like a series of strongly en- 
tangled pairs next to each other (with weaker links between) 
and this is shown for the = 6 case in step 1 of Fig. |8ja). 
When the system evolves, the state of the system takes the 
form of step 2 in Fig. |8ta) and after a while it goes to the 
form of step 3 in Fig. ISa). To explain this curious effect, 
without losing the generality, we consider the isotropic AFM. 
Clearly Fig. ISb), where a singlet between 0' and an odd site 
breaks 3 strong bonds, is energetically not favored in course of 
a unitary dynamics starting as step 1 of Fig. [8ja). So despite 
a finite (but small) overlap between the state shown in Fig. 
Htb) and those in Fig. |8la) this state does not emerge through 
the dynamics in the sense that its overlap with the state |V'(0) 
never become higher than a certain value. Quantitatively, all 
moments of Hamiltonian are conserved [90 during the evolu- 
tion: V71, = (?A(OI-ff"l^(0> = ( ^(0)|ij'"|V'(0)) , so 
energy (E = {H)) and its variance (rj = ^ {H'^) — {H)'^) are 
constant during the evolution. It means that only states with 
energy expectation E for which E — ri<E<E + 77, such as 
in Fig. [HI a) can contribute in evolution, while those as in Fig. 
HJb) cannot play a role. Note also that this curious phenom- 
ena, when recast in terms of two-time correlations, states that 
< crQ(0)cr|(<) > should be less than —1/3 only for the even 
sites j. Thus is potential physical systems where such dynam- 
ical correlations is measurable, the hopping mode of transfer 
should be testable. 



XI. POTENTIAL PHYSICAL REALIZATIONS 

We now mention some systems in which our results can be 
potentially tested, though there is some way to go for some of 
these systems, as local addressing of the spin to be suddenly 
coupled to the chain may be required. Recently there has been 
extensive interest in finite spin chains such as fabricated AFM 
nano-chains lfl2ll . and especially even-odd effects in such sys- 
tems lEsll . This can be one potential system where recently de- 
veloped sensitive magnetometers 112911 can perhaps be used for 
verifying the correlations and hence entanglement. Perhaps 
an STM tip encoding the spin to be transmitted can be brought 
close to one end of a finite array. Finite chains of dop ed ful- 
lerines in nanotubes (such as AFM Sc@C82 lEj) is the 
other alternative for developing this idea. Spins in such sys- 
tems have already been measured, and perhaps local electri- 
cal gates can give local control to couple in the input qubit 
ll30ll . Optical superlattices with atoms can realize an ensem- 
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ble of finite spin chains jsill as well as the switching on their 
interactions 1331 ■ Barrier heights at regular intervals may be 
raised to create arrays of small lattice segments (cells) of sizes 

2 and Net with the repeating pattern 2, Nch, 2, Nch, The 

O'O singlet and the finite chain ground states can be created in 
the cells of sizes 2 and Nch respectively as ground states (in 
fact, the former has already been accomplished ||3lll ). Next, 
again through global methods, the barriers between the 2 site 
cells and the Nch site cells to their right have to be lowered 
(simultaneously the barrier between the two sites of the cell 
of length 2 has to be raised), so as to form superlattices with 
cells of size Nch + 2 each. The subsequent dynamics will then 
exactly be as we have predicted and can potentially be ver- 
ified through global time of flight correlation measurements 
llsill . One can use ion traps where small spin systems are be- 
ing realized Issll . as well as implementing spin chains with 
trapped electrons llsill where initializing individual spins and 
controlling the interaction at one end are both simple. NMR 
is another fruitful avenue for testing communication through 
spin chains ifssll . 

XII. SUMMARY 

We have studied the transmission of both classical and 
quantum information through the all phases of the XX Z 
Hamiltonian. This quantifies the ability of each phase for 
information transmission. We found that in the absence of 
noise and thermal fluctuation isotropic Heisenberg Hamilto- 
nian (A = 1) is the best point of the phase diagram for in- 



formation transmission, both in terms of its amount, as well 
as its speed. The speed of propagation of the information, de- 
spite our finite open-ended case, fits strikingly well with the 
spin wave velocities known from continuum limit field theo- 
retic studies of the XXZ spin chain. When decoherence and 
thermal fluctuations are taken to account, the best point of the 
phase diagram moves to the Neel phase, which due to a faster 
evolution, is less sensitive to these sources of noise. Fur- 
thermore, we showed that the transmission through an even 
chain is characterized by the Pauli channel which has bene- 
fits in terms of immediate applicability of entanglement dis- 
tillation. We also studied the transmission of classical in- 
formation through this channel. Optimal states for single- 
shot classical capacity were identified and we realized that 
even when system is so noisy, such that quantum information 
is completely destroyed, some classical information can be 
transferred. Studying the entanglement propagation through 
the chain showed that entanglement skips odd numbered sites 
and manifests as a curious behavior of two time correlation 
functions during the non-equilibrium dynamics. It remains 
an open problem to explain well the mysterious behavior of 
the dynamics which entanglement suddenly drops around the 
point A = -0.5. 
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